AIAA 2002-1516 


BUCKLING BEHAVIOR OF COMPRESSION-LOADED COMPOSITE 
CYLINDERICAL SHELLS WITH REINFORCED CUTOUTS 

Mark W. Hilburger* and James H. Starnes, Jr.^ 

NASA Langley Research Center 
Hampton , Virginia 23681 -00 1 


Abstract 

Results from a numerical study of the response of 
thin-wall compression-loaded quasi-isotropic laminat- 
ed composite cylindrical shells with reinforced and un- 
reinforced square cutouts are presented. The effects of 
cutout reinforcement orthotropy , size, and thickness on 
the nonlinear response of the shells are described. A 
high-fidelity nonlinear analysis procedure has been 
used to predict the nonlinear response of the shells. The 
analysis procedure includes a nonlinear static analysis 
that predicts stable response characteristics of the shells 
and a nonlinear transient analysis that predicts unstable 
dynamic buckling response characteristics. The results 
illustrate how a compression-loaded shell with an unre- 
inforced cutout can exhibit a complex nonlinear re- 
sponse. In particular, a local buckling response occurs 
in the shell near the cutout and is caused by a complex 
nonlinear coupling between local shell- wall deforma- 
tions and in-plane destabilizing compression stresses 
near the cutout. In general, the addition of reinforce- 
ment around a cutout in a compression-loaded shell can 
retard or eliminate the local buckling response near the 
cutout and increase the buckling load of the shell, as ex- 
pected. However, results are presented that show how 
certain reinforcement configurations can actually cause 
an unexpected increase in the magnitude of local defor- 
mations and stresses in the shell and cause a reduction 
in the buckling load. Specific cases are presented that 
suggest that the orthotropy, thickness, and size of a cut- 
out reinforcement in a shell can be tailored to achieve 
improved response characteristics. 

Introduction 

Thin-walled cylindrical shells with cutouts are 
found in many aerospace structural applications. Ad- 
vanced material systems are being used to produce 
stronger, lighter- weight aerospace shell structures, and 
improved analysis and design methods appropriate for 
these advanced material systems are needed. The high 
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strength-to-weight and high stiffness-to-weight ratios 
of advanced composite materials offer significant 
weight reduction potential for aerospace structures. In 
addition, the use of advanced composite materials al- 
lows the designer to tailor the stiffness properties of the 
structure to obtain a structurally efficient design. Many 
aerospace shell structures have cutouts or openings that 
serve as doors, windows, or access ports and these cut- 
outs or openings often require some type of reinforcing 
structure to control local structural deformations and 
stresses near the cutout. In addition, these structures 
may experience compression loads during operation, 
and thus their buckling response characteristics must 
also be understood and accurately predicted in order to 
determine safe operating conditions and effective de- 
signs for these structures. 

Many studies have been conducted which show 
that a cutout in an isotropic shell structure can have a 
significant effect on the response of the shell. In partic- 
ular, results indicate that a cutout in a shell structure 
causes a local response to occur near the cutout when 
the shell is subjected to load. This local response can 
consist of large out-of-plane deformations, large-mag- 
nitude rapidly varying stresses near the cutout (e.g., 
Refs. 1-2). If the load is a compressive load, the cutout 
can cause a local buckling response to occur in the shell 
at applied loads lower than the general instability load 
of the corresponding shell without a cutout (e.g., Refs. 
3-5). For some cases, this local buckling response re- 
sults in a stable postbuckling response localized near 
the cutout and additional load can be applied to the shell 
before it exhibits global collapse. For other cases, the 
local response in the shell causes a disturbance in the 
shell of sufficient magnitude to cause global collapse to 
occur immediately after the local instability occurs. 

Significantly fewer studies have been conducted 
on the response of compression-loaded composite cy- 
lindrical shells with cutouts. Recent numerical and ex- 
perimental studies of the response of compression- 
loaded composite cylindrical shells with unreinforced 
rectangular cutouts have been presented by Hilburger, 

et al. 6, 7 In particular, their results have shown that lo- 
calized regions of biaxial in-plane compression stresses 
form in the shell near the cutout, and these regions of bi- 
axial stresses couple with the out-of-plane deforma- 
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tions of the shell wall causing an unstable local 
buckling response to occur near the cutout. In addition, 
their results indicate that laminate orthotropy can have 
a significant effect on the nonlinear response and the 
buckling behavior of a shell with a cutout, as expected. 
For example, numerical results from Ref. 6 indicate that 
the normalized buckling load* of a [±45/0/90] s quasi- 
isotropic shell with a 1 .0-in. by 1 .0-in. square cutout is 
approximately 33 % lower than the normalized buckling 
load of a corresponding [±45/0 2 ] s orthotropic shell with 
the same sized cutout. The large reduction in the nor- 
malized buckling load of the quasi-isotropic shell, as 
compared to the orthotropic shell, is caused by larger 
magnitude radial deformations and destabilizing in- 
plane compression stresses that develop near the cutout 
in the shell at a lower applied load level. These local- 
ized deformations and in-plane stresses cause the local 
buckling of the shell to initiate at a lower applied load 
level. These numerical results were verified with sev- 
eral experiments and the results indicated good agree- 
ment. 

Similarly, very little information is available on 
the response of compression-loaded curved shells with 
reinforced cutouts, and the few results that do exist are 
limited to isotropic shells (e.g.. Refs. 8-10). For exam- 
ple, Almroth and Holmes 8 presented results from a nu- 
merical and experimental study on the response of 
compression-loaded cylindrical shells with reinforced 
and unreinforced rectangular cutouts. Their results 
showed that the arrangement of the cutout reinforce- 
ment, that is, whether the reinforcement is positioned 
along an axially aligned free-edge or around all the edg- 
es of the cutout, can have a significant effect on the 
buckling response of the shell. Experimental results 

from Toda 10 indicate that, for certain sized circular cut- 
outs in a compression-loaded cylindrical shell, the 
bucking load increased with an increase in the size of 
the reinforcement around the cutout. However, further 
reinforcement of the cutout beyond a certain amount 
did not always produce higher buckling loads, and in 
some cases further reinforcement caused a decrease in 
the buckling load of the shell. 

A simplified design-level analysis procedure and 
empirical data are often used in the design of shells with 
reinforcements. The traditional method for the prelim- 
inary design of a reinforced cutout in a thin-walled shell 
structure is based on the linear analysis of a flat plate 


The normalized buckling load is defined as the 
buckling load of the shell with a cutout predicted 
from a nonlinear analysis normalized with respect to 
the buckling load of the corresponding shell without 
a cutout predicted from a linear bifurcation analysis. 


with a square cutout (e.g.. Ref. 1 1) and the effects of 
shell curvature and other modeling approximations are 
taken into account by applying empirical correction 
factors. However, this design method may not be suffi- 
cient for the design of large cutouts in curved shell 
structures because of the lack of accuracy in the analy- 
sis and the complex nonlinear behavior exhibited by 
these types of structures. In addition, this design meth- 
od can result in overly conservative designs or uncon- 
servative designs for these structures if the empirical 
correction factors are not applicable to the design of in- 
terest. For example, current design practice does not in- 
clude information or design procedures for the design 
of reinforced cutouts in composite shell structures. 

A review of the results presented in the literature 
indicates that the response of a compression-loaded cy- 
lindrical shell with an unreinforced cutout is, for the 
most part, understood. In contrast, the effects of cutout 
reinforcement on the buckling behavior of compres- 
sion-loaded composite cylindrical shells, is not well un- 
derstood. The objective of this paper is to identify the 
typical nonlinear response characteristics of a compres- 
sion-loaded, thin-walled, quasi-isotropic, laminated, 
cylindrical shell with a square cutout and to illustrate 
the effects of several cutout reinforcement configura- 
tions on the response. This thin-walled cylindrical shell 
configuration was chosen as a generic example of a typ- 
ical aerospace shell structure with a cutout subjected to 
a destabilizing load. Numerical results that show the ef- 
fects of 18 different reinforcement configurations on 
the response of these shell structures are presented. The 
various cutout reinforcement configurations considered 
were used to study the effects of reinforcement orthot- 
ropy , thickness, and size on the response of the shell. In 
addition, the results are used to illustrate the feasibility 
of structurally tailoring cutout reinforcement to control 
shell-wall deformations and stresse concentrations near 
a cutout in a compression-loaded shell and to increase 
the load carrying capacity of the shell. A modem high- 
fidelity nonlinear analysis procedure is used in the 
study (e.g.. Ref. 12) that offers the opportunity to pro- 
vide insight into the effects of various cutout reinforce- 
ment concepts on the response of compression-loaded 
shell structures. This procedure also offers the oppor- 
tunity to improve some of the engineering approxima- 
tions and methods that are used in the design of 
reinforced cutouts in shell structures. The high-fidelity 
nonlinear analysis procedure used in this study has been 
successfully applied to the analysis of other similar 
compression-loaded shells with cutouts and the predict- 
ed results have been verified with selected experiments, 
e.g., see Ref. 7 . First, results illustrating the response of 
a compression-loaded cylinder with an unreiforced cut- 
out are presented. Then, results illustrating the effects 
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of selected cutout reinforcement on the response of the 
shell are presented. Lastly, response trends are identi- 
fied and discussed. The results include load-shortening 
response curves; prebuckling, buckling and postbuck- 
ling deformation patterns; and prebuckling, buckling, 
and postbuckling stress contours. 

Finite-Element Models and Analyses 
Nonlinear Analysis Procedure 

The shells considered in this study were analyzed 
with the STAGS (STructural Analysis of General 

Shells) nonlinear shell analysis code. 13 STAGS is a fi- 
nite-element code designed for the static and dynamic 
analysis of general shells, and includes the effects of 
geometric and material nonlinearities in the analysis. 
The STAGS code uses both the modified and full New- 
ton methods for its nonlinear solution algorithms, and 
accounts for large rotations in a shell wall by using a co- 
rotational algorithm implemented at the element level. 

The Riks pseudo arc-length path-following method 14 is 
used to continue a solution past the limit points of a 
nonlinear response. With this strategy, the incremental- 
ly applied loading parameter is replaced by an arc- 
length along the solution path, which is then used as the 
independent loading parameter. The arc-length incre- 
ments are automatically adjusted by the program as a 
function of the solution behavior. A transient analysis 
option is included in STAGS that uses an implicit nu- 
merical time-integration method developed by Park 15 
and proportional structural damping. Mass and stiff- 
ness damping factors used in the transient analysis are 
defined as cc=2jcv£, and P=t/2jw, respectively , where £ 
denotes the fraction of critical damping and is assumed 
to equal 0.15 for these shells. The dynamic buckling 
response of the shell is assumed to be dominated by a 
particular frequency of vibration, v, and is defined as 
the critical frequency of vibration determined from a 
linear vibration analysis in STAGS. Additional infor- 
mation on the transient analysis in STAGS can be found 
in Ref. 16. 

The prebuckling, buckling and postbuckling re- 
sponses of the shells were determined by using the fol- 
lowing analysis procedure. The quasi-static responses 
were determined by using the geometrically nonlinear 
quasi-static analysis capability in STAGS. The Riks 
pseudo arc-length path-following method was used to 
compute the response of the shell up to the buckling 
point. The unstable buckling response of the shell was 
predicted by using the nonlinear transient analysis op- 
tion of the code. The transient analysis was initiated 
from an unstable equilibrium state close to the buckling 
point by incrementing the end displacement by a small 
amount, approximately 1%. An initial time step of 10" 8 


seconds was used in the analysis and is automatically 
adjusted by the program as a function of the solution be- 
havior. The transient analysis was continued until the 
kinetic energy in the shell dissipated to a negligible lev- 
el, which indicates that the transient response has atten- 
uated. Once the transient analysis had attenuated to a 
near-steady-state solution, the load relaxation option of 
the code was used to establish a stable static equilibri- 
um state. Conventional linear bifurcation buckling 
analysis results were also determined with STAGS for 
comparison with the nonlinear response results. 

Finite-Element Models 

Nineteen shells were analyzed in the present study 
and include one shell with an unreinforced cutout and 
18 shells with reinforced cutouts; these shells are re- 
ferred to herein as shells Cl through C 19. A list of shell 
identification codes and corresponding reinforcement 
configurations is given in Table 1. The shells consid- 
ered in the present study are modeled as geometrically 
perfect 8-ply laminated quasi-isotropic circular cylin- 
ders with a square cutout and subjected to a uniform ax- 
ial end-shortening. The shells have a length of 16.0 in., 
a radius of 8.0 in., and shell- wall thickness of 0.04 in. 
that give a shell-radius-to-thickness ratio of 200. The 
cutout is a 1 .0 in. by 1 .0 in. square-shaped cutout with 
0.05-in-radius reentrant comers. The shell wall was 
modeled as a [±45/0/90] s quasi-isotropic laminate, in 
which each lamina ply had a thickness of 0.005 in. and 
a fiber volume fraction of 0.62. The material properties 
of a lamina ply are as follows: longitudinal compression 
modulus E 1 = 18.5 Msi, transverse modulus E 2 = 1.64 
Msi, in-plane shear modulus G 12 = 0.87 Msi, and major 
Poisson’s ratio v i2 =0.30. Shells C2 through C19 in- 
clude various forms of reinforcement around the cutout. 
The cutout reinforcement consists of additional square- 
shaped lamina plies added to the shell-wall laminate at 
the shell-wall mid-surface and are aligned concentrical- 
ly with respect to the square-shaped cutout in the shell. 
Eighteen different reinforcement configurations were 
considered in which the reinforcement thickness, size, 
and ply-orientation were varied. Specifically, three dif- 
ferent reinforcement thickness of0.005 in., 0.01 in., and 
0.02 in., corresponding to one, two, and four-ply rein- 
forcement thicknesses, respectively, were considered. 
Two reinforcement ply orientations of 0° and 90° (a 0° 

ply and a 90° ply correspond to lamina plies with fibers 
aligned along the length of the shell and around its cir- 
cumference, respectively) were considered. Three rein- 
forcement sizes: a 2.4-in. by 2.4-in. square 
reinforcement, a 4.4-in. by 4.4-in. square reinforce- 
ment, and a 8.0-in. by 8.0-in. square reinforcement. 

A typical finite-element model of a specimen is 
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shown in Fig. 1 . The standard 410 quadrilateral ele- 
ment from the STAGS element library was used in the 
models. Typically, the elements of the finite-element 
mesh are approximately 0.2-in. by 0.2-in. square, but 
the mesh is refined near the cutout. Each element pos- 
sesses four integration points, which are distributed to 
provide a modeling resolution of at least 0.1 -in. Results 
from a mesh convergence study indicated that this de- 
gree of mesh refinement mesh was required to accurate- 
ly model the bending boundary layer deformation 
response near the ends of the shell, the prebuckling 
short wavelength deformation response and rapidly 
varying stress distributions exhibited by the shells near 
the cutout, and the short wavelength deformation re- 
sponse and rapidly varying stress distributions exhibit- 
ed by the shells during the transient buckling process. 
Clamped end conditions were used in the finite-element 
models in which the circumferential and radial dis- 
placements v and w are set equal to zero in end regions 
of the shell as indicated in Fig. 1 . A typical finite-ele- 
ment model contained approximately 120,000 degrees 
of freedom. 

Results and Discussion 

Numerically predicted results for selected com- 
pression-loaded quasi-isotropic laminated cylindrical 
shells with reinforced and unreinforced cutouts are pre- 
sented in this section. The results were obtained from 
finite-element models of geometrically perfect shells 
subjected to a uniform axial end-shortening. These re- 
sults are presented to illustrate the overall behavior of a 
compression-loaded graphite-epoxy shell with a cutout 
and the effects of cutout reinforcement on the response. 
First, results illustrating the nonlinear response of a 
compression-loaded geometrically perfect quasi-isotro- 
pic cylindrical shell with an unreinforced square- 
shaped cutout are presented. Then, results illustrating 
the predicted response of selected compression-loaded 
cylindrical shells with reinforced cutouts are presented 
and compared. The results include load-shortening re- 
sponse curves; prebuckling, buckling and postbuckling 
deformation patterns; and prebuckling, buckling, and 
postbuckling stress contours. 

Cylindrical Shell with an Unreinforced Cutout 

Results for a quasi-isotropic shell with a 1 .0-in by 
1 .0-in unreinforced square-shaped cutout, referred to 
herein as shell Cl , are presented in this section to illus- 
trate typical response characteristics of a compression- 
loaded shell with a square cutout. The numerically pre- 
dicted load-shortening response curve for shell Cl is 
shown in Fig. 2a. The quasi-static response is indicated 
by a solid line and the unstable, transient buckling re- 
sponses are indicated by the dashed lines in the figure. 


The axial load P is normalized with respect to the bifur- 
cation buckling load predicted from a linear eigenvalue 
analysis of the corresponding shell without a cutout, P cr 
= 42,590 lbf, and the end-shortening displacement A is 
normalized with respect to the shell length, L. The shell 
exhibits a complex nonlinear response as indicated by 
the numerous stable and unstable segments of the load- 
shortening response curve. The shell has a linear preb- 
uckling response and a distinct buckling point labeled A 
at a normalized load level of P/P cr = 0.51. The buckling 
response is characterized by a local, unstable transient 
buckling response in the shell and includes the forma- 
tion of large magnitude radial deformations and rapidly 
varying stresses in the shell wall near the cutout. The 
shell exhibits a relatively small reduction in axial load, 
approximately 6%, during buckling and is indicated by 
the dashed line-segment A-B of the load-shortening re- 
sponse curve. The figure also presents a second solu- 
tion path from point A to point B . This path 
corresponds to the quasi-static solution predicted by the 
arc-length path-following method. The results indicate 
that the quasi-static analysis and the transient analysis 
predict the same local postbuckling state at point B. 
During the prebuckling response, regions of destabiliz- 
ing in-plane biaxial compression stresses form near the 
cutout and couple with the radial deformations near the 
cutout to cause the unstable local buckling response in 
the shell. A stable postbuckling response is exhibited 
by the shell as indicated by segment B-C of the load- 
shortening response. As loading continues in the post- 
buckling region of the response, the magnitude of the 
shell- wall radial deformations near the cutout increases . 
In addition, the results indicate that the slope of the 
postbuckling curve decreases as loading continues in 
the postbuckling range, indicating a reduction in effec- 
tive axial stiffness. This decrease in effective axial 
stiffness is caused by increasingly large radial deforma- 
tions that develop in the shell that causes a significant 
redistribution of load away from the cutout to reduce 
the effective, load-carrying cross-section of the shell. 
The global collapse of the shell occurs at point C and 
the unstable collapse response is represented in Fig. 2b 
by the dashed line segment C-F. The collapse response 
is characterized by a significant reduction in load and 
the development of the general instability mode. The 
transient nature of the unstable buckling responses in 
the shell is indicated by the load-time history of the col- 
lapse response of the shell shown in Fig. 2b. The load 
history curve exhibits a sudden reduction in axial load 
supported by the shell from P/P cr = 0 .5 1 , at time t = 0 .0 
seconds, to P/P cr = 0.32, at time t = 0.006 seconds into 
the transient response. In addition, the kinetic energy in 
the shell increases rapidly during the initial collapse re- 
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sponse and then dissipates as time progresses until the 
shell reaches a stable postbuckling equilibrium config- 
uration. The shell obtains a stable postbuckling config- 
uration after approximately 0.015 second have elapsed 
and the normalized postbuckling load value for the shell 
is equal to 0.33. 

Deformation patterns corresponding to selected 
points during the compression response are shown in 
Figs. 3a-3f, and deformation patterns with superim- 
posed axial and circumferential stress contours corre- 
sponding to selected points during the compression 
response are shown in Figs. 4a-4c (the magnitude of the 
deformation responses in these figures are exaggerated 
for clarity). The deformation response incipient to local 
buckling, (associated with point A in Fig. 2a), shown in 
Fig. 3a, is characterized by two ellipse-shaped inward 
buckles near the cutout with the semi-major axis of the 
buckles aligned in a helical or skew direction. The 
magnitude of the shell-wall radial deformations vary 
between -0.65 and +0.40 times the shell-wall thickness, 
where positive- and negative- valued displacements cor- 
respond to outward and inward radial displacements, 
respectively. The skewing of the deformation pattern is 
attributed to the presence of laminate anisotropy in the 
form of coupling between the bending and twisting of 
the shell wall. The axial and circumferential membrane 
stress contours indicate significant stress concentra- 
tions near the cutout that rapidly decay to far-field val- 
ues away from the cutout, as shown in Fig. 4a 
(magnified views of the stresses near the cutout are 
shown for clarity). The maximum value of the axial 
compression stress resultant occurs at the comers of the 
cutout (-2023 lb/in.) and is equal to 4.7 times the corre- 
sponding uniform far-field value. Similarly, large mag- 
nitude circumferential stress concentrations occur near 
the cutout with stress resultant values that range be- 
tween -695 lb/in. and +528 lb/in; the maximum circum- 
ferential tension stress resultant occurs at the mid-point 
of the upper and lower free-edges of the cutout and the 
maximum circumferential compression stress resultant 
occurs at the cutout comers. In addition, regions of lo- 
calized circumferential compression stress resultant de- 
velop approximately 0.5 inches above and below the 
cutout with maximum stress resultant levels of approx- 
imately -110 lb/in. The local buckling response of the 
shell is initiated by a nonlinear coupling between local- 
ized destabilizing compressive axial and circumferen- 
tial stress resultants and the radial deformations that 
occur in the shell near the cutout (e.g., Fig. 4a). The lo- 
cations where local buckling initiates are marked by X 
symbols in Fig. 4a. The initial postbuckling deforma- 
tion pattern associated with point B in Fig. 2a is shown 
in Fig. 3b and is characterized by large ellipse-shaped 
buckles on either side of the cutout. The magnitude of 


the shell-wall radial deformations have increased sig- 
nificantly and vary between -4.5 and +1.5 times the 
shell-wall thickness. As loading continues in the post- 
buckling range, the ellipse-shaped buckles in the shell 
wall increase in size and the buckle pattern rotates 
clockwise around the cutout, as shown in Fig. 3c. The 
magnitude of the shell-wall radial deformations in- 
creases to between -8.65 and +4.2 times the shell- wall 
thickness and a significant redistribution of the load in 
the shell occurs near the buckles, as shown in Fig. 4b. 
More specifically, the axial stress resultant distribution 
indicates a significant unloading of the shell near the in- 
ward buckles. This load redistribution coincides with 
the reduction in postbuckling stiffness indicated in Fig. 
2a and causes a reduction in the effective, load-carrying 
cross-section of the shell. Moreover, a significant in- 
crease in the magnitude of the circumferential compres- 
sion stress resultants occurs in the buckles. At the onset 
of global collapse, the shell-wall deformations change 
in a very short period of time from the stable, local de- 
formation pattern shown in Fig. 3c to a unstable tran- 
sient deformation pattern shown in Fig. 3d. Moreover, 
additional buckles form in the shell wall around the cir- 
cumference of the shell. As the buckling process con- 
tinues, there is a significant reduction in the axial load 
supported by the shell, and even more buckles develop 
in the shell wall, as indicated in Fig. 3e. After approx- 
imately 0.015 seconds have elapsed, the kinetic energy 
in the shell has dissipated to a negligible level, and the 
shell has deformed into a stable post-collapse general 
instability mode-shape that consists of eight buckles 
distributed around the circumference of the shell as 
shown in Fig. 3f. The corresponding post-collapse ax- 
ial and circumferential stress resultant distribution pat- 
terns are shown in Fig. 4c and indicate significant stress 
resultant redistribution throughout the entire cylinder. 
In particular, the results indicate that the majority of the 
applied axial compression load is now supported by the 
ridges between the adjacent inward buckles in the shell 
wall, and alternating bands of circumferential tension 
and compression stress resultants are present in the 
shell. 

Effects of Cutout Reinforcement 

Results for eighteen compression-loaded quasi- 
isotropic shells with 1 .0-in by 1 .0-in square reinforced 
cutouts are presented in this section to identify the ef- 
fects of selected cutout reinforcement configurations on 
the response of the shells (identified herein as shells 
C2-C19). The results for these shells are compared to 
results for the corresponding shell with an unreinforced 
cutout, shell Cl, presented in the previous section. 

The buckling results for the compression-loaded 
shell with an unreinforced cutout, shell Cl , presented in 
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the previous section, identified several features of the 
behavior that are associated with the local buckling re- 
sponse of the shell. Most significantly, the nonlinear 
interaction between localized radial deformations and 
destabilizing biaxial stresses near the cutout cause the 
local buckling response of the shell to occur. This fun- 
damental behavioral feature suggests that it may be pos- 
sible to identify a cutout reinforcement configuration 
that retards the onset of the local deformations and de- 
stabilizing stress in the shell near the cutout and can re- 
tard or eliminate the onset of the local buckling 
response in the shell. To this end, a numerical paramet- 
ric study was conducted to identify the effects of select- 
ed cutout reinforcement configurations on the response 
of a compression-loaded shell with a cutout. In partic- 
ular, cutout reinforcement size, thickness, and orthotro- 
py were studied. All of the reinforcements considered 
herein consist of square-shaped concentrically aligned 
lamina plies added to the shell-wall laminate at the 
shell-wall mid-surface. Thus, three reinforcement sizes 
were investigated to identify the effects of reinforce- 
ment size in the deformation response and stress distri- 
bution in the shell. The three square-shaped 
reinforcement sizes include 2.4-in. by 2.4-in, 4.4-in. by 
4. 4- in, and 8.0-in by 8.0-in square reinforcement (re- 
ferred to herein as 2.4-in-square reinforcement, 4.4-in- 
square reinforcement, and 8.0-in-square reinforcement, 
respectively). Three reinforcement thicknesses were 
studied and include reinforcement thickness equal to 
0.005 in., 0.01 in., and 0.02 in., corresponding to one-, 
two-, and four-ply reinforcements, respectively. Two 
reinforcement ply orientations of 0° and 90° (a 0°-ply 

and a 90°-ply correspond to lamina plies with fibers 
aligned along the length of the shell and around its cir- 
cumference, respectively) were investigated to study 
the effects of reinforcement orthotropy on the response. 
A list of shell identification codes and corresponding 
reinforcement configurations is given in Table 1 . 

2.4-in. Square Reinforcement . Results from 
shells C2 through C7 are presented in this section. The 
cutout reinforcements for shells C2 through C4 and 

shells C5 through C7 consist of 0° and 90° 2.4-in. 
square-shaped lamina plies, respectively. Shells C2 
and C5 have 0.005-in-thick reinforcement, shells C3 
and C6 have 0 .0 1 -in-thick reinforcement, and shells C4 
and C7 have 0.02-in- thick reinforcement. Numerically 
predicted load-shortening response curves for shells C2 
through C4 and shells C5 through Cl, are presented in 
Figs. 5a and 5b, respectively, and the results for the cor- 
responding unreinforced shell, Cl , are shown for com- 
parison. The quasi-static responses are indicated by the 
solid lines and the unstable, transient buckling respons- 


es are indicated by the dashed lines in the figures. The 
axial load P is normalized with respect to the bifurca- 
tion buckling load of the corresponding shell without a 
cutout predicted from a linear eigenvalue analysis P cr = 
42,590 lbf, and the end-shortening displacement A is 
normalized with respect to the shell length L. The local 
buckling points are marked with filled circles, and glo- 
bal collapse points are indicated by an X. Curves show- 
ing the unstable transient responses associated with 
local buckling (e.g., line segment A-B shown in Fig. 2a 
for shell Cl) for shells C 1 , C2, and C5 are omitted from 
Figs. 5a and 5b for clarity. The results show that the ad- 
dition of reinforcement around the cutout can have a 
significant effect on the overall character of the load- 
shortening response of the shell. In particular, the re- 
sults indicate that the local buckling loads can increase 
significantly as the thickness of the reinforcement is in- 
creased. More specifically, the normalized buckling 
loads for shells C2, C3, and C4 are equal to 0.57, 0.64, 
and 0.78, respectively, as compared to the buckling 
load of 0.51 for the corresponding unreinforced shell 
C 1 . Similarly, the normalized buckling loads for shells 
C5 , C6 , and C7 are equal to 0 .59 , 0 .67 , and 0 .77 , respec- 
tively . However, the results indicate that some of the 
shells with thicker reinforcement, shells C3, C4, C6, 
and C7, no longer exhibit stable post-local-buckling re- 
sponses. For these shells, the local buckling response in 
these shell causes enough of a disturbance in the shell 
to cause the global collapse of the shell. In contrast, 
shells C2, and C4 exhibit stable post-local-buckling re- 
sponse and can support additional load before global 
collapse occurs. Upon global collapse, each shell ob- 
tains a stable post-collapse configuration and exhibits a 
significant reduction in axial load-carrying capability. 
The normalized post-collapse load levels for shell C2- 
C7 range from 0.32 to 0.35. For the most part, the re- 
sults presented in Figs. 5a and 5b indicate that the cut- 
out reinforcement orientation, that is 0°-ply orientation 
versus 90°-ply orientation, has a small effect (at most a 
5% difference) on the buckling load of the shell and a 
negligable effect on the effective axial stiffness of the 
shell for this size reinforcement. 

Displacement contour patterns just before buck- 
ling for shells C2-C4 are shown in Figs. 6a-6c, respec- 
tively . The radial displacements w are normalized with 
respect to the shell -wall thickness t = 0.04 in. The outer 
boundaries of the cutout reinforcement in each shell is 
marked in the figures with solid lines. The results indi- 
cate that the character of the deformation response in 
the shell incipient to buckling changes significantly as 
the thickness of the cutout reinforcement increases. In 
particular, the displacement contour patterns exhibited 
by shells C2 and C3 in Figs. 6a and 6b, respectively, in- 
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dicate similar local deformation response as that exhib- 
ited by shell Cl which is characterized by two ellipse- 
shaped inward buckles near the cutout. The semi-major 
axis of the buckles are aligned in a helical direction and 
the localized displacement gradients rapidly attenuate 
away from the cutout. However, the magnitudes of the 
shell-wall deformations decrease as the thickness of the 
reinforcement increases; that is, the magnitude of the 
normalized displacements range between -0.51 to 0.3 
and -0.33 to 0.29 for shells C2 and C3, respectively, as 
compared to -0 .65 to 0 .4 for shell C 1 . These results cor- 
respond to as much as 21% and 49% reductions in the 
magnitude of the local deformations near the cutout for 
shells C2 and C3 , respectively , as compared to shell C 1 . 
The deformation response for shell C4 just before buck- 
ling, shown in Fig. 6c, exhibits a further reduction in the 
displacement gradients in the shell wall near the cutout, 
with an increase in the thickness of the reinforcement, 
as compared to shells C2 and C3 . More specifically, the 
magnitude of the shell wall displacements range be- 
tween -0.06 and 0.33 and correspond to a 91% reduc- 
tion in the magnitude of the maximum inward 
displacement near the cutout as compared to shell Cl . 
In addition, the maximum inward shell-wall displace- 
ment in shell C4 occurs just outside the reinforced re- 
gion of the shell, and this response is unlike that of 
shells C1-C3 in which the maximum inward shell-wall 
displacement occurs at the free-edges of the cutout. 
Axial and circumferential traces illustrating the effects 
of the cutout reinforcement on the deformation re- 
sponse in shells Cl through C4 are compared in Figs. 7a 
and 7b, respectively. The regions of the shell that the 
trace data were taken from is marked by the dashed ax- 
ial and circumferential lines superimposed on the dis- 
placement contour plot shown in Fig. 6a. The results 
show that the magnitude of the shell-wall displace- 
ments decrease near the cutout as the reinforcement 
thickness is increased, as expected. Furthermore, the 
results indicate that the character of the attenuation of 
the local displacements in shell C4 are somewhat differ- 
ent from those exhibited by shells Cl through C3. In 
particular, shell C4 exhibits additional waves in the dis- 
placement pattern along the length and around the 
circumference of the shell. 

The character of the axial and circumferential 
stress resultant distributions just before buckling for 
shells C2 through C4 are similar to that of shell Cl 
shown in Fig. 4a. In particular, the shells exhibit regions 
of in-plane biaxial compression stress resultants near 
the cutout, however, the magnitudes of the stresses can 
be significantly different. More specifically, the mag- 
nitudes of the axial and circumferential stress resultant 
components in these regions of biaxial compression just 
before buckling are as follows: N x = -580 lb/in and N e 


= -109 lb/in for shell C2, N x = -615 lb/in and N 0 = -96 
lb/in for shell C3, and N x = -703 lb/in and N e = -87.5 lb/ 
in for shell C4, as compared to N x = -568 lb/in, and N e 
= -110 Ib/in for shell Cl . For comparison purposes, it 
is useful to define an axial stress concentration factor 
and a biaxial-stress ratio. The axial-stress concentra- 
tion factor is defined as the local axial stress resultant 
value in the region of biaxial stress, N x , normalized 
with respect to the far-field value of axial stress result- 
ant N x °= P/2jtR, where P is the axial load applied to the 
shell and R is the shell radius. The biaxial-stress ratio 
is defined as the ratio of local circumferential stress to 
axial stress in the region of biaxial compression, that is 
N 0 /N x . The definitions of stress concentration factor 
and biaxial stress ratio are used herein unless otherwise 
specified. The stresses in the region of biaxial compres- 
sion correspond to axial stress concentration factors 
equal to 1.20, 1.1 3, and 1.06, for shells C2,C3,andC4, 
respectively, as compared to 1 .3 1 for shell Cl , and bi- 
axial stress ratios of 0.19, 0.16, and 0.12 for shells C2, 
C3, and C4, respectively, as compared to 0.19 for shell 
Cl . These results indicate that the axial stress concen- 
tration factor and the biaxial-stress ratio decrease as the 
thickness of the cutout reinforcement increases. 

The results indicate that the overall character of 
the local buckling response in shells C2-C4 is similar to 
that of shell Cl . In particular, the local buckling re- 
sponses in shells C2-C4 initiate near the cutouts and are 
caused by a nonlinear coupling between the localized 
regions of destabilizing in-plane biaxial compression 
stress resultants and the local shell-wall deformations. 
In addition, the locations at which buckling initiates in 
shells C2-C4 correspond to points within the reinforced 
region of the shells. However, the results also indicate 
that the character of the post-local-buckling response 
and the transient collapse response can be greatly af- 
fected by the reinforcement. For example, shell C2 ex- 
hibits a stable post-local-buckling response and global 
collapse response similar to that of shell Cl . However, 
when the reinforcement thickness is increased, the 
shells no longer exhibit a stable post-local-buckling re- 
sponse; rather, the shell exhibits a global collapse, e.g., 
shells C3 and C4, and the shells collapse into a general 
instability mode-shape. The general instability mode 
shape for shell C3 is similar to that exhibited by shell 
Cl (see Fig. 3f), however, shell C4 buckles into a mode 
shape that consists of two half waves along the length 
of the shell and eight full waves around the circumfer- 
ence. Upon local buckling and global collapse, shells 
C2-C4 exhibit significant stress redistribution through- 
out the entire shell, similar to that exhibited by shell Cl 
as shown in Figs. 4b, and 4c. 

Displacement contour patterns just before buck- 
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ling for shells C5-C7 (corresponding shells with 90°- 
ply cutout reinforcement) are shown in Figs. 8a-8c, re- 
spectively. The radial-displacement contour patterns 
for shells C5 through C7 indicate large-magnitude local 
deformations that rapidly attenuate away from the cut- 
out, and are similar to the contours exhibited by shell 
Cl (see Fig. 3a). The magnitudes of the shell-wall de- 
formations decrease as the thickness of the reinforce- 
ment increases, and the magnitude of the normalized 
displacements range from -0.58 to 0.35 for shell C5, - 
0.54 to 0.38 for shell C6, and -0.46 to 0.42 for shell C7, 
as compared to -0.65 to 0.4 for shell Cl . These results 
correspond to a reduction in the magnitude of the local 
deformations by as much as 29% as the thickness of the 
reinforcement increases, as compared to shell Cl. 
However, the reduction in the magnitude of the local 
displacements for these shells is significantly less than 
that exhibited by the corresponding shells with 0°-ply 
reinforcement which exhibited a maximum reduction 
of approximately 91%. In addition, axial and circum- 
ferential displacement traces illustrating the effects of 
the cutout reinforcement on the deformation response 
of shells C5 through C7 are compared in Figs. 9a and 
9b, respectively. These results indicate that the thick- 
ness of the cutout reinforcement has a small effect on 
the magnitude and distribution of the shell wall defor- 
mations near the cutout. This behavior is in contrast to 
the significant variation in the deformation responses 

exhibited by the corresponding shells with 0°-ply rein- 
forcement (shells C2 through C4) as shown in Figs. 7a 
and 7b. 

The character of the axial and circumferential 
stress resultant distributions just before buckling for 
shells C5 through C7 are, for the most part, similar to 
that of shell Cl shown in Fig. 4a. However, the magni- 
tudes of the stresses are significantly different. In par- 
ticular, the magnitudes of the axial and circumferential 
stress resultant components in these regions of biaxial 
compression near the cutout just before buckling are as 
follows: N x = -620 lb/in and N 0 = -103 lb/in for shell 
C5, N x = -662 Ib/in and N 0 = -94 Ib/in for shell C6, and 
N x = -759 lb/in and N 0 = -89 lb/in for shell C7, as com- 
pared to N x = -568 lb/in, and N 0 = -1 10 lb/in for shell 
Cl. These results correspond to axial-stress concentra- 
tion factors equal to 1 .25, 1 .19, and 1 .16, for shells C5, 
C6, and C7, respectively, as compared to 1 .31 for shell 
Cl and biaxial-stress ratios of 0.17, 0.14, and 0.12 for 
shells C5, C6, and C7, respectively, as compared to 
0.19 for shell Cl. These results indicate that the axial- 
stress concentration factor and the biaxial-stress ratio 
decrease as the thickness of the cutout reinforcement in- 
creases. This trend is similar to the response trend ex- 


hibited by shells C2-C4. In addition, the results for 
these shells indicate that the magnitude of the axial- 
stress concentration factor and the biaxial stress ratio 
are, on average, 6.2% greater and 8.7% lower, respec- 
tively, as compared to the corresponding shells with 0°- 
ply reinforcement, shells C2-C4. 

The results indicate that the overall character of 
the local buckling response in shells C5-C7 is similar to 
that of shell Cl . In addition, the locations at which the 
buckling response initiates in the shells correspond to 
points within the reinforcement regions. However, the 
results also indicate that the character of the post-local- 
buckling response and the transient collapse response 
can be affected by reinforcement thickness. For exam- 
ple, shell C5 exhibits a similar stable post-local-buck- 
ling and global collapse responses as that of shell Cl 
and C2. However, as the reinforcement thickness is in- 
creased, the shell no longer exhibits a stable post-local- 
buckling response; rather, the local buckling response 
in the shell causes the global collapse of the shell to oc- 
cur (e.g., shells C6 and C7). In addition, shells C5 and 
C6 exhibit post-collapse general instability mode 
shapes similar to that of shell Cl . However, as the 
thickness of the reinforcement is increased, the shell ex- 
hibits a general instability mode shape that consists of 
two half waves along the length of the shell and eight 
full waves around the circumference of the shell. This 
mode shape is similar to the post-collapse mode shape 
exhibited by shell C4. Overall, these response trends 
are similar to those exhibited by the corresponding 
shells with 0°-ply cutout reinforcements, shells C2-C4. 
Upon local buckling and collapse, shells C5-C7 exhibit 
significant stress redistribution throughout the entire 
shell, similar to that exhibited by shell Cl that is shown 
in Figs. 4b, and 4c. 

4.4-in. Square Reinforcement . Results for shells 
C8 through C13 are presented in this section. The cut- 
out reinforcements for shells C8 through C 10 consist of 

4.4-in square-shaped 0° lamina plies and the cutout re- 
inforcements for shells Cl 1 through C13 consist of 4.4- 

in square-shaped 90° lamina plies. Shells C8 and Cl 1 
have 0.005-in- thick reinforcement, shells C9 and C12 
have 0.01-in-thick reinforcement, and shells CIO and 
C13 have 0.02-in-thick reinforcement. Numerically 
predicted load-shortening response curves for shells C8 
through CIO and shells Cl 1 through Cl 3, are presented 
in Figs. 10a and 10b, respectively. The quasi-static re- 
sponses are indicated by the solid lines and the unstable, 
transient buckling responses are indicated by the 
dashed lines in the figure. The axial load P is normal- 
ized with respect to the bifurcation buckling load of the 
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corresponding shell without a cutout P cr = 42,590 lbf, 
and the end-shortening A is normalized with respect to 
the shell length L. The local buckling points are 
marked with filled circles, and global collapse points 
are indicate by an X. 

The results in Figs. 10a and 10b also show that the 
addition of reinforcement around the cutout can have a 
significant effect on the overall character of the re- 
sponse of the shell, and indicate similar trends to those 
shown by shells C2-C7 in Figs. 5a and 5b. In particular, 
the local buckling loads increase significantly as the 
thickness of the reinforcement increases. Specifically, 
the normalized buckling loads for shells C8, C9, and 
CIO are equal to 0.62, 0.73, and 0.84, respectively, as 
compared to the buckling load of 0.51 for the corre- 
sponding unreinforced shell Cl . Similarly, the normal- 
ized buckling loads for shells Cl 1 , C12, and C13 are 
equal to 0.64, 0.76, and 0.92, respectively. The results 
also indicate that the shells with 90°-ply reinforcements 
considered here can exhibit significantly higher buck- 
ling loads than the corresponding shells with 0°-ply re- 
inforcements, that is, the buckling loads for shells Cl 1 , 
C12, and C13 are 3.2, 4.1, and 9.5% higher than for 
shells C8, C9, and CIO, respectively. In addition, the 
results indicate that the buckling loads of the shells in- 
crease as the reinforcement size increases. More spe- 
cifically, the buckling loads for shells C8, C9, and CIO 
are 10, 18, and 12% higher than the buckling loads for 
the corresponding shells with 2.4-in square-shaped 0°- 
ply reinforcements , shells C2 , C3 , and C4 , respectively . 
Similarly, the buckling loads for shells C14-C16 are 9, 
18, and 29% higher than the corresponding buckling 
loads for shell with 2.4-in-square 90°-ply reinforce- 
ments C5, C6, and C7, respectively. Furthermore, the 
results indicate that some of the shells with thicker re- 
inforcement, C9 through Cl 3, no longer exhibit stable 
post-local-buckling responses. Rather, the local buck- 
ling response in each shell causes a disturbance of suf- 
ficient magnitude to cause the global collapse of the 
shell. This response characteristic is contrary to that of 
shell C8, which exhibits a stable post-local-buckling re- 
sponse before global collapse, similar to that exhibited 
by shell Cl (see Fig. 2a). Upon global collapse, each 
shell obtains a stable post-collapse configuration, ac- 
companied by a significant reduction in axial load-car- 
rying capability. The normalized post-collapse load 
levels range from 0.31 to 0.35 and these post-collapse 
load levels are similar to those exhibited by the corre- 
sponding shells with smaller reinforcement; that is, 
shells C2-C7 (see Figs. 5a and 5b). 

Displacement contour patterns just before buck- 
ling for shells C8-C10 are shown in Figs. 1 la-1 lc, re- 
spectively. The radial displacements w are normalized 


with respect to the shell-wall thickness t =0.04 in. The 
outer boundary of the cutout reinforcement in each shell 
is marked in the figures with solid lines. In general, the 
results indicate that the character of the deformation 
patterns just before local buckling can change signifi- 
cantly as the thickness of the reinforcement increases. 
In addition, the magnitude of the local deformations 
near the cutout decreases as the thickness of the rein- 
forcement increases. In particular, the displacement 
contour patterns exhibited by shell C8, shown in Fig. 

1 la, exhibit a large-magnitude local deformation re- 
sponse that rapidly attenuates away from the cutout in a 
manner that is similar to that of shell Cl . The magni- 
tude of the normalized shell wall deformations range 
between -0 .47 to 0.35 for shell C8 , as compared to -0.65 
to 0.4 for shell Cl , and correspond to as much as 28% 
reduction in the magnitude of the local deformations 
near the cutout, as compared to shell Cl . However, the 
character and the magnitude of the local deformations 
change significantly as the thickness of the reinforce- 
ment increases. In particular, the magnitude and char- 
acter of the displacement gradients near the cutout in 
shell C9 just before buckling are similar to the displace- 
ment gradients exhibited by shell C8. However, the 
magnitude of the deformations away from the cutout 
and the attenuation lengths of these deformations in- 
crease as the thickness of the reinforcement increases, 
as shown in Figs. 1 lb, 12a, and 12b. Furthermore, the 
results indicate that the displacement gradients near the 
cutout are reduced significantly as the thickness of the 
cutout reinforcement is increased further and exhibit 
normalized displacement magnitudes that range from 
0.37 to 0.57; e.g., see shell CIO shown in Fig. 1 lc. In 
addition, the large-magnitude displacements away from 
the cutout, and the attenuation length of these displace- 
ments increase, and the shell exhibits normalized dis- 
placement values as large as -0.40, as shown in Figs. 
11c, 12a, and 12b. Overall, these results indicate that 
the local deformations near the cutout decrease with an 
increase in the thickness of the cutout reinforcement. 
However, the shells also exhibit increasing large-mag- 
nitude deformations away from the cutout and the atten- 
uation length of these deformations increases as the 
thickness of the reinforcement is increased. This re- 
sponse trend is similar to that exhibited by the corre- 
sponding shells with the smaller 2.4-in-square 
reinforcement, however, the magnitudes of the defor- 
mations away from the cutout are as much as four times 
larger. Additionally, the attenuation lengths are ap- 
proximately twice as long in shells C8-C10, as those ex- 
hibited by the corresponding shells with the smaller 
2.4-in-square reinforcement. 

The character of the local axial and circumferen- 
tial stress resultant distributions just before buckling for 
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shells C8 through CIO are, largely, similar to those pre- 
dicted for shell Cl shown in Fig. 4a. However, the 
magnitudes of the stress resultants can vary significant- 
ly as the thickness of the cutout reinforcement increas- 
es. In particular, the magnitudes of the axial and 
circumferential stress resultant components in these re- 
gions of biaxial compression just before buckling are 
N x = -719 lb/in and N 0 = -129 lb/in for shell C8, and N x 
= -865 lb/in and N 0 = -126 lb/in for shell C9 (N x = -568 
lb/in, and N e = -1 10 lb/in for shell Cl). These results 
correspond to axial-stress concentration factors equal to 
1 .36 and 1 .40 for shells C8 and C9, respectively, as 
compared to 1 .31 for shell Cl and biaxial-stress ratios 
of 0.18 and 0.15 for shells C8 and C9, respectively, and 
0.19 for shell Cl . In addition, the magnitude of the ax- 
ial-stress concentration factor and the biaxial stress ra- 
tio for shells C8-C10 are, 22, 5, and 5.7% greater on 
average than those exhibited by the corresponding 
shells with the smaller sized cutout reinforcement, C2- 
C4, respectively. In contrast, the axial and circumfer- 
ential stress resultant distributions change significantly 
as the thickness of the reinforcement increases further. 
Specifically, shell CIO exhibits localized regions of in- 
plane biaxial compression stress resultants in the bend- 
ing boundary layer region of the shell and the magni- 
tudes of the stress are N x = -985 lb/in and N 0 = -58 lb/ 
in. These stress resultants correspond to an axial-stress 
concentration factor of 1 .3 8 , and a biaxial-stress ratio of 
0.06. 

The deformation patterns and stresses indicate 
that the initiation of the local buckling response of 
shells C8 and C9 is similar to that of shell Cl . Howev- 
er, the results also indicate that the character of the post- 
local-buckling response and the transient collapse re- 
sponse can be greatly affected by the cutout reinforce- 
ment. In particular, shell C8 exhibits a similar stable 
post-local buckling response and global collapse re- 
sponse as that exhibited by shell Cl as shown in Figs. 
3c-3f. However, as the thickness of the reinforcement 
increases the shell no longer exhibits a stable post-lo- 
cal-buckling response. Instead, the local response in 
shell causes the overall global collapse to occur (e.g., 
shell C9). This response trend for shells C8 and C9 is 
similar to that exhibited by the corresponding shells 
with smaller reinforcement, shells C2 and C3. As the 
thickness of the reinforcement is increased further, the 
shell no longer exhibits a local buckling response near 
the cutout. For this case, a local buckling response ini- 
tiates in the bending boundary layer region where large- 
magnitude inward deformations occur, e.g., see shell 
CIO, and Fig. 11c). These deformations couple with lo- 
calized destabilizing in-plane biaxial compression 
stress resultants in the bending boundary layer region of 


the shell. This localized buckling response in the bend- 
ing boundary region of the shell causes the global col- 
lapse of the shell. This response characteristic is unlike 
that of the corresponding shell with the smaller, 2.4-in- 
square cutout reinforcement (shell C4) which exhibits 
local buckling near the cutout. This result illustrates 
how the reinforcement size can have a significant effect 
on the character of the buckling response of the shell. 
In particular, the cutout reinforcement has suppressed 
the local buckling response near the cutout and the shell 
now exhibits a different localized response in the bend- 
ing boundary region that causes the buckling of the 
shell. Upon local buckling and collapse, shells C8-C10 
exhibit significant stress redistribution throughout the 
entire shell in a manner similar to that of shell Cl that 
is shown in Figs. 4b, and 4c. 

The deformation response just before buckling for 
shells Cl 1-C13 (corresponding shells with 90°-ply cut- 
out reinforcement) are shown in Figs. 13a-13c, respec- 
tively and also illustrate how the thickness of the 
reinforcement can have a significant effect on the defor- 
mation response. In particular, the radial-displacement 
contour patterns of shells C 1 1 and C 1 2 show similar lo- 
cal large-magnitude deformations that rapidly attenuate 
away from the cutout and are similar to those exhibited 
by shell C 1 . However, as the thickness of the reinforce- 
ment increases , the shell exhibits significant differences 
in the deformation pattern (e.g., shell C13), as com- 
pared to shells Cl 1 and Cl 2. Specifically, large-mag- 
nitude radial deformations develop in a large portion of 
the shell away from the cutout and include two large in- 
ward buckles that develop in the bending boundary lay- 
er region of the shell. In addition, the magnitudes of the 
shell wall deformations increase as the thickness of the 
reinforcement increases from 0.005 inches to 0.01 inch- 
es; that is, the magnitudes of the shell-wall deforma- 
tions range between -0.67 to 0.39 for shell Cl 1 and 
between -0.96 to 0.68 for shell C12, as compared to - 
0.65 to 0.4 for shell Cl . These results correspond to an 
increase in the magnitude of the displacements by as 
much as 70% in shells Cl 1 and C12. However, upon 
further increase of the reinforcement thickness to 0.02 
inches, the magnitude of the normalized shell- wall de- 
formations decrease significantly to between -0.37 to 
0 .45 for shell C 1 3 , and correspond to a 43 % decrease in 
the magnitude of the deformations, as compared to shell 
Cl . Axial and circumferential traces illustrating the ef- 
fects of the cutout reinforcement on the deformation re- 
sponse of shells C 1 1 through C 1 3 are compared in Figs . 
14a and 14b, respectively. The results show that signif- 
icant changes in the character of the deformation re- 
sponse of the shells occur as the thickness of the 
reinforcement is changed. In particular, shell C 13 ex- 
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hibits a significant increase in the attenuation length of 
the local deformation response as compared to shells 
Cl 1 and Cl 2. These results indicate that, for the most 
part, the character and magnitudes of the deformation 
responses for shells with 4.4-in-square 90°-ply rein- 
forcement (see shells Cl 1 through C13 shown in Figs. 
1 la-1 lc) can be significantly different from the corre- 
sponding shells with 0°-ply reinforcement (see shells 
C8 through CIO in Figs. 13a-13c). In particular, the re- 
inforcement in shells C8-C10 reduces the local defor- 
mations near the cutout just before buckling. In 
contrast, shells Cl 1-C13 can exhibit increases in the lo- 
cal deformations with an increase in the reinforcement 
thickness. 

The character of the local axial and circumferen- 
tial stress resultant distributions near the cutout just be- 
fore buckling for shells Cl 1 through C13 are, for the 
most part, similar to that of shell Cl shown in Fig. 4a. 
However, the magnitudes of the stress resultants can be 
significantly different. In particular, the magnitudes of 
the axial and circumferential stress resultant compo- 
nents in these regions of biaxial compression near the 
cutout just before buckling are as follows: N x = -681 lb/ 
in and N e = -126 Ib/in for shell Cl 1 , and N x = -843 lb/ 
in and N 0 = -173 lb/in for shell C12, as compared to N x 
= -568 lb/in, and N 0 = -1 10 lb/in for shell Cl . These re- 
sults correspond to axial-stress concentration factors 
equal to 1.26 and 1.30 (1.31 for shell Cl) and biaxial- 
stress ratios of 0.19 and 0.21 , for shells Cl 1 and C12, 
respectively (0.19 for shell Cl). These results indicate 
that the axial-stress concentration factor and the biaxial 
stress ratio increase as the thickness of the cutout rein- 
forcement increases. The magnitude of the axial-stress 
concentration factor and the biaxial-stress ratio for 
these shells are, on average, approximately 7.2% lower 
and 21% greater, respectively, as compared to the cor- 
responding shells with 0°-ply reinforcement (shells C8- 
C9). This trend exhibited by the shells with 4.4-in- 
square reinforcements is contrary to that of the shells 
with smaller 2.4-in-square cutout reinforcements, in 
which the shells with 0°-ply reinforcement exhibited 
larger-magnitude axial-stress concentration factors and 
lower biaxial-stress ratios than that of the correspond- 
ing shells with 90° reinforcement. In addition, the re- 
sults indicate that the magnitude of the stress 
concentration factors and the biaxial-stress ratios are, 
on average, 0.8 and 39.9% greater than those exhibited 
by the corresponding shells with the smaller reinforce- 
ment, shells C5-C7. In contrast, the axial and circum- 
ferential stress resultant distributions change 
significantly as the thickness of the reinforcement in- 
creases further to 0.02 inches. Specifically, shell C13 


exhibits localized regions of in-plane biaxial compres- 
sion stress resultants in the bending boundary layer re- 
gion of the shell and the magnitudes of the stress 
resultants are N x = -838 lb/in, and N e = -128 lb/in, 
which correspond to an axial-stress concentration factor 
of 1 .08, and a biaxial-stress ratio of 0.15. This response 
characteristic is similar to that of shell CIO. The results 
also indicate that shells C11-C13 exhibit similar post- 
local-buckling and collapse response trends as those ex- 
hibited by the corresponding shells with 0°-ply rein- 
forcement, shell C8-C10. 


8.0-in. Square Reinforcement . Results for shells 
C14 through C19 are presented in this section. The cut- 
out reinforcements for shells C14 through Cl 6 consist 

of 8.0-in by 8.0-in square-shaped 0° lamina plies and 
the cutout reinforcements for shell C17 through C19 
consist of 8.0-in by 8.0-in square-shaped 90° lamina 
plies. Shells C 14 and C 17 have 0.005-in-thick rein- 
forcement, shells C15 and C18 have 0.01-in-thick rein- 
forcement, and shells C16 and C19 have 0.02-in-thick 
reinforcement. Predicted load-shortening response 
curves for shells C14 through C16 and shells C17 
through C19, are presented in Figs. 15a and 15b, re- 
spectively. The quasi-static responses are indicated by 
solid lines and the unstable, transient buckling respons- 
es are indicated by dashed lines in the figure. The axial 
load P is again normalized with respect to the bifurca- 
tion buckling load of the corresponding shell without a 
cutout obtained from a linear eigenvalue analysis P cr = 
42,590 lbf, and the end-shortening displacement A is 
normalized with respect to the shell length L. The local 
buckling points are marked with filled circles, and glo- 
bal collapse points are indicated by an X. These results 
show that the addition of reinforcement around the cut- 
out can have a significant effect on the overall character 
of the response of the shell, and exhibit, for the most 
part, similar response trends to those shown by shells 
C2-C7 in Figs. 5a and 5b, and shells C8-C13 in Figs. 
10a and 10b. In particular, the local buckling loads of 
the shells increase as the thickness of the reinforcement 
increases; that is, the normalized buckling loads for 
shells C14, C15, and C16 are equal to 0.61, 0.72, and 
0.74, respectively, as compared to the normalized buck- 
ling load of 0.51 for the corresponding unreinforced 
shell Cl . However, the results indicate that an increase 
in the size of the cutout reinforcement in shells 04- 
06 causes a 1 .6, 1 .4, and 12% reduction in the buck- 
ling loads as compared to the corresponding buckling 
loads for shells C8-C10 (shells with 4.4-in-square 0°- 
ply reinforcements), respectively. Furthermore, the re- 
sult indicate that the effective axial compression stiff- 
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ness of the shell increases as the thickness of the 
reinforcement increases, as shown in Fig. 15a. This re- 
sult is in contrast to the results from shells C2-C4 and 
C8-C10 which exhibit negligible changes in the effec- 
tive axial stiffness of the shells with the addition of the 
cutout reinforcements. The normalized buckling loads 
for shells C17,C18,andC19 equal 0 .66 , 0 .83 , and 0 .93 , 
respectively, and correspond to a 3.1, 9.2, and 1% in- 
crease in the buckling load over the corresponding 
buckling loads for shells Cl 1 , C12, and C13 (shells 
with 4.4-in-square 90°-ply reinforcements), respective- 
ly. The results indicate that the shells with 90°-ply re- 
inforcements considered herein exhibit significantly 
higher buckling loads than the corresponding shells 

with 0°-ply reinforcements; that is, the buckling loads 
for shells C17, C18, and C19 are 8.2, 15.3, and 25.7% 
higher than the buckling loads of shells C14, C15, and 
Cl 6, respectively. This response trend is similar to that 
exhibited by shells C8-C 13, however, the differences in 
the buckling loads are significantly greater for the shells 
with the 8.0-in-square reinforcement. Furthermore, the 
results indicate that some shells with thicker cutout re- 
inforcement no longer exhibit stable post-local-buck- 
ling responses. Rather, the local buckling response in 
each shell causes a disturbance in the shell of sufficient 
magnitude to cause the global collapse of the shell, e.g., 
shells C15 through C19. In contrast, shell C14 which 
exhibits a stable post-local-buckling response before 
global collapse similar to that exhibited by shell Cl (see 
Fig. 2a). Upon global collapse, each shell obtains a sta- 
ble post-collapse configuration and exhibits a signifi- 
cant reduction in axial load carrying capability. The 
normalized post-collapse load levels for shells Cl 4- 
C19 range from 0.34 to 0.39 and are as much as 26% 
greater than the corresponding post-collapse load-lev- 
els of shells C2-C7 and C8-C13. 

Displacement contour patterns just before buck- 
ling for shells C14-C16 are shown in Figs. 16a-16c, re- 
spectively. The radial displacements w are normalized 
with respect to the shell- wall thickness t = 0.04 in. The 
outer boundary of the cutout reinforcement in each shell 
is marked in the figures with solid lines. These results 
also indicate that the character of the deformation re- 
sponse in the shell just before buckling changes signif- 
icantly as the thickness of the reinforcement increases. 
For example, the displacement contour patterns exhib- 
ited by shells C14 and C15 shown in Figs. 16a and 16b, 
respectively, exhibit similar local large-magnitude de- 
formations that rapidly attenuate away from the cutout 
and are similar to the response exhibited by shell Cl . In 
addition, the magnitudes of the normalized shell-wall 
deformations increase as the thickness of the reinforce- 
ment increases and range between -0.54 to 0.48 for 


shell C14 and between -0.98 to 0.62 for shell C15, and 
correspond to as much as a 20% and 50.8% increase in 
the magnitude of the local deformations near the cutout, 
respectively, as compared to shell Cl . In contrast, the 
character of deformation response just before buckling 
for shell Cl 6, changes significantly with an increase in 
the thickness of the reinforcement as shown in Fig . 1 6c , 
as compared to shells C14 and C15. In particular, the 
local deformation response near the cutout is character- 
ized by a large ellipse-shaped inward buckle and the 
semi-major axis of the buckle is aligned with the longi- 
tudinal axis of the shell. In addition, large-magnitude 
outward radial deformations occur on either side of the 
cutout. Furthermore, the displacement gradients near 
the cutout in shell C16 are reduced significantly with 
normalized displacement magnitudes that range from - 
0.30 to 0.47, e.g., see Figs. 17a and 17b, and correspond 
to a 27% reduction in the local displacement gradient. 
Generally, the results indicate that the magnitude of the 
shell-wall displacements away from the cutout and the 
attenuation length of these shell wall displacements in- 
crease as the thickness of the reinforcement increases, 
as shown in Figs. 17a and 17b, respectively. Further- 
more, the magnitude of the deformations and the atten- 
uation length of these deformations are significantly 
larger than those exhibited by the corresponding shells 
with 2.4-in- and 4.4-in-square cutout reinforcement. 

The character of the local axial and circumferen- 
tial stress resultant distributions near the cutout just be- 
fore buckling for shells C14 through C16 are, for the 
most part, similar to those of shell Cl (shown in Fig. 
4a). However, the magnitudes of the stress resultants 
can vary significantly as the thickness of the reinforce- 
ment increases. The magnitudes of the axial and cir- 
cumferential stress resultant components in the regions 
of bi-axial compression just before buckling are N x = - 
719 lb/in and N 0 = -121 lb/in for shell C14, and N x = - 
940 lb/in and N 0 = -227 Ib/in for shell C15. The results 
correspond to axial-stress concentration factors equal to 
1.39, 1.55, for shells C14,C15, respectively, (1.31 for 
shell Cl) and biaxial-stress ratios of 0.17, 0.24, for 
shells C 14, Cl 5, respectively, (0.19 for shell Cl). The 
stress concentration factors for these shells are, on av- 
erage, 26.2 and 6.5% greater than those exhibited by the 
corresponding shells with 2.4-in- and 4.4-in-square re- 
inforcements, respectively. Similarly, the results indi- 
cate that biaxial-stress ratios for these shells are, on 
average, 20.0 and 27.3% greater than those exhibited by 
the corresponding shells with 2.4-in- and 4.4-in-square 
reinforcements, respectively. In contrast, shell C16 de- 
veloped regions of biaxial stress resultants in the bend- 
ing boundary region of the shell with stress resultant 
values of N x = -867 lb/in, and N 0 = -163 lb/in for shell 
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Cl 6, and correspond to an axial-stress concentration 
factor of 1 .38 and a biaxial-stress ratio of 0.16. This re- 
sponse characteristic is similar to that exhibited by the 
corresponding shell with 4.4-in-square reinforcement, 
shell CIO. The results indicate that the initiation of the 
local buckling response of shells C14 and C15 is similar 
to that of shell Cl, however, the results also indicate 
that the character of the post-local-buckling response 
and the transient collapse response can be effected by 
the reinforcement. For example, shell C14 exhibits a 
similar stable post-local-buckling response and global 
collapse response as that of shell Cl . As the thickness 
of the reinforcement increases, the shell no longer ex- 
hibits a stable post-local-buckling response. Rather, the 
local buckling response in the shell causes the overall 
collapse of the shell (e.g., shell C15). As the thickness 
of the reinforcement is increased further, the shell no 
longer exhibits a local buckling response near the cut- 
out, rather, a local buckling response initiates in the 
bending boundary layer region of the shell similar to 
that exhibited by shell CIO. This localized buckling re- 
sponse in the bending boundary layer region of the shell 
initiates the global collapse of the shell. This response 
trend associated with the post-local-buckling and glo- 
bal collapse response is similar to that exhibited by the 
corresponding shells with 4.4-in-square reinforcement. 

Displacement contour patterns just before buck- 
ling for shells C 17-09 (corresponding shells with 90°- 
ply cutout reinforcement) are shown in Figs. 18a-18c, 
respectively, and indicate that the cutout reinforcement 
can have a significant effect on the deformation re- 
sponse in the shell. In particular, the radial displace- 
ment contour patterns exhibited by shells C17 and Cl 8 
indicate local large-magnitude deformations similar to 
those for the corresponding 0°-ply reinforced shells 
04 and 05 and the unreinforced shell Cl . The mag- 
nitudes of the shell- wall deformations range from -0.8 
to 0.41 for shell 07 and -0.71 to 0.45 for shell 08, and 
indicate that the magnitude of the local deformations in- 
crease by as much as 23 and 9.2% in shells 07, 08, 
respectively, as compared to shell Cl . However, as the 
thickness of the reinforcement increases further, the 
shell exhibits a significant reduction in the local defor- 
mations (e.g., shell 09). In addition, large-magnitude 
radial deformations develop in a large portion of the 
shell away from the cutout and include four large in- 
ward ellipse-shaped buckles that develop in the bending 
boundary layer region of the shell and is a similar defor- 
mation response to that of shell 06 (e.g., Fig. 16c). 
The magnitude of the shell-wall deformations ranges 
from -0.37 to 0.36 and correspond to a 43% increase in 
the deformations away from the cutout. Axial and cir- 
cumferential traces illustrating the effects of the cutout 


reinforcement on the deformation response of shells 
07 through 09 are compared in Figs. 19a and 19b, re- 
spectively. The results show the significant changes in 
the character of the deformation response of the shells 
as the thickness of the reinforcement is changed. In par- 
ticular, shell 09 exhibits a slight increase in the mag- 
nitude of the shell-wall deformations away from the 
cutout and an increase in the attenuation length of the 
local deformation response around the circumference of 
the shell as compared to shells 07 and 08. These re- 
sults indicate that, for the most part, the character of the 
deformation responses for shells with 8.0-in-square 

90°-ply reinforcement (e.g., see shells 07 through 09 
Figs. 18a- 18c), is somewhat similar to that of the corre- 
sponding shells with 0°-ply reinforcement (e.g., see 
shells 04 through 06 in Figs. 16a- 16c). However, 
the magnitudes and attenuation lengths of the deforma- 
tions in shells 07-09 are generally smaller then those 
of shells 04-06. In addition, the magnitude of the 
deformations in shells with 8 .0-in-square reinforcement 
can be as much as 38% larger than those of the corre- 
sponding shells with 2.4-in-square reinforcement, and 
26% smaller than those exhibited by the corresponding 
shells with 4.4-in-square reinforcement. 

The character of the local axial and circumferen- 
tial stress resultant distributions just before buckling for 
shells 07 through 09 are, for the most part, similar to 
that of shell Cl (shown in Fig. 4a). However, the mag- 
nitudes of the stress resultants can be significantly dif- 
ferent as the thickness of the reinforcement increases. 

In particular, the axial and circumferential stress result- 
ant components in these regions of biaxial compression 
just before buckling are N x = -736 lb/in and N 0 = -156 
lb/in for shell 07, and N x = -888 Ib/in and N e = -145 
lb/in for shell 08. These results correspond to axial- 
stress concentration factors equal to 1.32 and 1.27 for 
shells 07 and 08, respectively, (1.31 for shell Cl) 
and biaxial-stress ratios of 0 .2 1 and 0 . 1 6, for shells C 1 7 
and 08, respectively (0.19 for shell Cl). These results 
indicate that the axial-stress concentration factor and 
the biaxial-stress ratio decrease as the thickness of the 
cutout reinforcement increases. In addition, the results 
indicate that the magnitude of the axial-stress concen- 
tration factors for these shells is, on average, 1 .9% low- 
er and the biaxial-stress ratio is, on average, 9.8% lower 
than those exhibited by the corresponding shells with 

0°-ply reinforcement. Furthermore, the results indicate 
that the stress concentration factors are, on average, 0.8 
and 1 .7% greater than the corresponding shells with 
2 ,4-in- , and 4 ,4-in-square reinforcements , respectively . 
Similarly, the results indicate that the biaxial-stress ra- 
tios are, on average, 30.6% greater and 10.5% less than 
the corresponding shells with 2.4-in-, and 4.4-in-square 
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reinforcements, respectively. However, as the thick- 
ness of the reinforcement increases further, the shell ex- 
hibits regions of biaxial compression stress resultants in 
the bending boundary layer region of the shell similar to 
that exhibited by shell C16 and the stress resultants are 
N x = -865 lb/in and Nq = -96 lb/in and correspond to a 
stress concentration factor of 1 .09 and a biaxial-stress 
ratio of 0. 1 1 . In general, the results indicate that shells 
C17-C19 exhibit similar post-local-buckling and col- 
lapse response trends as those exhibited by the corre- 
sponding shells with 0°-ply reinforcement, shell 04- 
06. 

Response Trends and Design Considerations 

The results presented herein illustrate how a com- 
pression-loaded composite shell with an unreinforced 
cutout can exhibit a complex nonlinear response, and 
that the character of the response can be significantly 
affected by the addition of reinforcement around the 
cutout. The buckling response in the shell is caused by 
a complex nonlinear coupling between local shell-wall 
deformations and in-plane destabilizing compression 
stresses near the cutout. The cutout reinforcement can 
retard the onset of buckling by retarding the onset of the 
local deformations and stress concentrations in the shell 
near the cutout. The results presented herein indicate 
that, for the most part, as the thickness and the size of 
the cutout reinforcement increases, the buckling load of 
the shell increases, as expected. In addition, the results 
show that reinforcement orientation can have a very 
pronounced effect on the shell response. In particular, 
shells with 90°-ply cutout reinforcement have higher 
buckling loads than the corresponding shells with 0° 
cutout reinforcement. Furthermore, shells with 0°-ply 
reinforcement exhibit smaller-magnitude deformations 
near the cutout than the corresponding shells with 90°- 

ply reinforcement. However, the shells with 0°-ply re- 
inforcement can exhibit larger-magnitude deformations 
away from the cutout, as compared to the correspond- 
ing shells with 90° reinforcement. In particular, as the 
size and thickness of the reinforcement increases, the 
magnitude and attenuation length of the deformations 
away from the cutout increases. These changes in char- 
acter of the deformation response can be important in 
understanding the formation of local bending gradients 
in thin-walled laminated shells. In particular, large- 
magnitude bending gradients in thin-walled laminated 
shells have been shown to cause interlaminar material 
failures (e.g., see Ref. 7). The results also indicate that, 
in general, shells with 1 or 2-ply reinforcement, i.e., 
0.005-in- and 0.0 1-in-thick reinforcement, exhibit local 
buckling near the cutout. The local buckling of the 


shell is caused by a complex nonlinear coupling be- 
tween local shell-wall deformations and in-plane desta- 
bilizing compression stresses near the cutout, and this 
response is similar to that exhibited by the correspond- 
ing shell with an unreinforced cutout, Cl . In contrast, 
some shells with 4-ply reinforcements no longer exhibit 
local buckling near the cutout. Rather, the shells exhib- 
it local buckling in the bending boundary layer region 
near the ends of the shell. For these cases, localized de- 
stabilizing in-plane biaxial stresess form in the bending 
boundary region of the shell and these stresses couple 
with the bending boundary layer deformations to cause 
local buckling to initiate in the bending boundary layer 
region of the shell. The results also indicate that shells 
with 1-ply reinforcement tend to exhibit stable post-lo- 
cal-buckling responses and can support additional load 
before global collapse. In contrast, the local buckling 
response of shells with thicker reinforcement (2-, and 4- 
ply-thick) can cause the overall global collapse of the 
shells. The results also indicate that shells with thinner 
reinforcements exhibit post-collapse general instability 
mode-shapes similar to that of the corresponding shell 
with an unreinforced cutout and are characterized by 
one half wave along the length of the shell and eight full 
waves around the circumference. In contrast, shells 
with thicker reinforcement can exhibit a post-collapse 
general instability mode-shapes that consists of two 
half waves along the length and eight full waves around 
the circumference. 

The results indicate that shells with reinforced cut- 
outs exhibit local stress distributions near the cutout 
similar to that exhibited by the corresponding shell with 
an unreinforced cutout. In particular, the shells exhibit 
significant axial and circumferential stress concentra- 
tions near the cutout and regions of in-plane biaxial 
compression stress. However, the magnitudes of the 
stress concentrations and the biaxial-stress ratio just be- 
fore buckling can vary significantly as the reinforce- 
ment orientation, size, and thickness vary. In particular, 
shells with 2.4-in-square cutout reinforcements exhibit 
axial-stress concentration factors that are, on average, 
11% lower than the corresponding shell with an unrein- 
forced cutout. Similarly, these shells exhibit biaxial- 
stress ratios that are, on average, 25% greater than that 
of the shell with an unreinforced cutout. In addition, 
there is at most a 5% difference in the stress concentra- 
tion factors and the biaxial-stress ratios exhibited by 
shells C2-C7. As the reinforcement size increase, the 
axial-stress concentration factors and the biaxial-stress 
ratios increase significantly for shells with 0°-ply rein- 
forcements and the stress concentration factor and the 
biaxial-stress ratio are, on average, 12.2% and 70.8% 
greater than that of the unreinforced shell Cl, respec- 
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tively. Similarly, shells with 90°-ply reinforcement ex- 
hibit biaxial-stress ratios that are, on average, 54% 
greater than that of the shell with an unreinforced cut- 
out. However, the axial-stress concentration factors re- 
main , on average ,1.1% lower than that exhibited by the 
unreinforced shell. These results suggest that shells 
with the 90°-ply reinforcement continue to achieve 
higher buckling loads as the reinforcement size increas- 
es, as compared to the corresponding shells with 0°-ply 
reinforcement, because the shells with 90°-ply rein- 
forcement do not generate increasingly larger magni- 
tudes of axial-stress concentrations as the 
reinforcement size increases. In contrast, shells with 

0°-ply reinforcement exhibit significant increases in the 
magnitude of the radial deformations and the axial- 
stress concentrations near the cutout at lower applied 
load levels. 

For the most part, the response trends for the shells 
presented herein are consistent with what one might ex- 
pect, however, the shells with 8.0-in-square 0°-ply cut- 
out reinforcement, C 14-06, exhibited some 
unexpected behavioral characteristics. In particular, 
these shells exhibited lower buckling loads than those 
exhibited by the corresponding shells with smaller 4.4- 
in-square reinforcements. Most notably, shell Cl 6, 
with a 4-ply-thick reinforcement, exhibited a normal- 
ized buckling load of 0.74 as compared to 0.84 for the 
corresponding shell with a 4.4-in square cutout rein- 
forcement, shell CIO. Furthermore, these shells exhibit 
substantially larger magnitudes of radial deformations 
and larger magnitudes of stress concentrations and bi- 
axial-stress ratios near the cutout than those exhibited 
by corresponding shells with smaller-size reinforce- 
ments. The results indicate that the load in the shell is 
redistributed into the region of the shell near the rein- 
forcement and causes the higher magnitude destabiliz- 
ing stresses and deformations to occur in the bending 
boundary layer region of the shell at lower applied load 
levels and ultimately cause the buckling of these shells 
to occur at lower applied load levels. In addition, these 
shells exhibit buckling loads that are significantly lower 
than the corresponding linear bifurcation buckling 
loads. In particular, the predicted buckling load for 
shell C16 is 12.6% lower than the corresponding bifur- 
cation buckling load. 

Concluding Remarks 

The results from a numerical study of the response 
of thin-wall compression-loaded quasi-isotropic lami- 
nated composite cylindrical shells with reinforced and 
unreinforced square cutouts have been presented. The 
results identify some of the effects of cutout-reinforce- 
ment orthotropy, size, and thickness on the nonlinear 


response of the shells. A high-fidelity nonlinear analy- 
sis procedure has been used to predict the nonlinear re- 
sponse of the shells. 

The results presented herein indicate that a com- 
pression-loaded shell with an unreinforced cutout can 
exhibit a complex nonlinear response that is not typical- 
ly found in the corresponding shell without a cutout. In 
particular, a local buckling response occurs in the shell 
that is caused by a localized nonlinear coupling be- 
tween local shell-wall deformations and in-plane desta- 
bilizing compression stresses near the cutout. In 
general, the addition of reinforcement around a cutout 
in a compression-loaded shell can have a significant ef- 
fect on the shell response, as expected. In particular, re- 
sults have been presented that indicate that the 
reinforcement can affect the local deformations and 
stresses near the cutout so as to retard or suppress the 
onset of local buckling in the shell near the cutout. For 
some cases, the local buckling response in the shell re- 
sults in a stable postbuckling response near the cutout 
and additional load can be applied to the shell before it 
undergoes a global collapse. For other cases, the local 
response in the shell causes a disturbance of sufficient 
magnitude in the shell to cause global collapse to occur 
immediately after the local instability occurs. For still 
other cases, the reinforcement suppresses the local 
buckling response near the cutout and causes buckling 
to occur in the bending boundary region of the shell. 
For the most part, the buckling load of the shell increas- 
es as the size and thickness of the cutout reinforcement 
increases. However, some results indicate that shells 
with a larger 0°-ply reinforcement can actually exhibit 
an increase in the magnitude of local deformations and 
stresses in the shell, and cause a reduction in the buck- 
ling load as compared to the corresponding shell with a 
smaller reinforcement. For these cases, the results indi- 
cate that this reinforcement configuration causes a re- 
distribution of the loads into the reinforced region of the 
shell; causes larger-magnitude displacements and 
stresses to occur at lower applied-load levels in the 
shell; hence, causes the buckling of the shell at lower 
applied load levels. The results indicate that shells with 

90°-ply cutout reinforcement tend exhibit higher buck- 
ling loads than the corresponding shells with 0°-ply re- 
inforcement. The results suggest that shells with the 
90°-ply reinforcement continue to achieve higher buck- 
ling loads because the shells with 90°-ply reinforce- 
ment exhibit smaller-magnitude axial-stress 

concentrations. In contrast, shells with 0°-ply rein- 
forcement exhibit larger-magnitude axial-stress con- 
centrations near the cutout. 

Overall, the results suggest that tailoring the 
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orthotropy , thickness and size of a cutout reinforcement 
in a compression-loaded shell can result in significant 
increases in the buckling load, and can reduce the local 
deformations and stresses near the cutout. However, 
the results also suggest that a nonlinear analysis may be 
necessary to accurately characterize the complex non- 
linear behavior exhibited by these structures. A modem 
high-fidelity nonlinear analysis procedure was used in 
the study and offers the opportunity to provide insight 
into the effects of various cutout reinforcement con- 
cepts on the response of compression-loaded shell 
structures. Moreover, results from a high-fidelity anal- 
ysis procedure can improve some of the engineering ap- 
proximations and methods that are used in the design of 
reinforced cutouts in shell structures. 
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Table 1. Shell Identification Codes and Reinforcement Configuration 


Shell Identification 
Code 

Reinforcement 
Size, 
in. x in. 

Reinforcement Ply 
Orientation, deg. 

Reinforcement 
thickness, in. 
(# of plies) 

Cl 

- 

- 

- 

C2 

2.4 x 2.4 

0 

.005 (1) 

C3 

2.4 x 2.4 

0 

.01 (2) 

C4 

2.4 x 2.4 

0 

.02(4) 

C5 

2.4 x 2.4 

90 

.005 (1) 

C6 

2.4 x 2.4 

90 

.01 (2) 

Cl 

2.4 x 2.4 

90 

.02 (4) 

C8 

4.4 x 4.4 

0 

.005 (1) 

C9 

4.4 x 4.4 

0 

.01 (2) 

CIO 

4.4 x 4.4 

0 

.02 (4) 

Cll 

4.4 x 4.4 

90 

.005 (1) 

C12 

4.4 x 4.4 

90 

.01 (2) 

C13 

4.4 x 4.4 

90 

.02 (4) 

C14 

8.0 x 8.0 

0 

.005 (1) 

C15 

8.0 x 8.0 

0 

.01 (2) 

C16 

8.0 x 8.0 

0 

.02 (4) 

C17 

8.0 x 8.0 

90 

.005 (1) 

C18 

8.0 x 8.0 

90 

.01 (2) 

C19 

8.0 x 8.0 

90 

.02 (4) 
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z, w 



region, v = w = 0 


Fig. 1 Typical finite-element model geometry and boundary conditions. 



Fig. 2 Numerically predicted nonlinear response of a compression-loaded cylindrical shell with an unrein- 
forced cutout Cl. 
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d) Deformation Pattern D e) Deformation Pattern E f) Deformation Pattern F 

Time - 0 .0022 seconds Time = 0 .0035 seconds Time = 0 .0 1 5 seconds 

P/P cr = 0.47 P/P cr = 0.39 P/P cr = 0.33 


Fig. 3 Numerically predicted deformation response of a compression-loaded cylindrical shell with an unrein' 
forced cutout Cl. The deformation patterns A-F correspond to the points labeled in Fig. 2. 
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magnified view magnified view 

a) Initial local buckling, corresponding to Point A in Fig. 2a, P/P cr = 0.51. 



b) Initial global collapse, corresponding to Point C in Fig. 2a, P/P cr = 0.56, time = 0.0 seconds. 



c) Post collapse, corresponding to Point F in Fig. 2a, P/P cr = 0.33, time = 0.015 seconds. 


Fig. 4 Numerically predicted deformation response and stress resultant contours at selected points in the non- 
linear response of a compression-loaded cylindrical shell with an unreinforced cutout Cl. 
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Fig. 5 Effects of cutout-reinforcement thickness and reinforcement fiber orientation on the predicted load- 
shortening response of a geometrically perfect compression-loaded cylindrical shell with a cutout. (2.4-in 
square-shaped reinforcement) 



mewi i mm-M n 

-0.51 w/t 0.30 


a) 0.005-in-thick reinforcement 



-0.33 w/t 0.29 

b) 0.01 -in-thick reinforcement 



ED J 1 1 IvIMTTl 

-0.06 w/t 0.33 

c) 0.02-in-thick reinforcement 


Fig. 6 Effects of cutout-reinforcement thickness on the predicted initial-buckling deformations of a geometri- 
cally perfect compression-loaded cylindrical shell with a cutout. The reinforcement is a 2.4-in. square rein- 
forcement with 0°-ply orientation (the outer boundary of the reinforcement region is marked in the figures). 
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Cl, no reinforcement 
C2, 0.005-in-thick reinforcement 
C3, 0.01 -in-thick reinforcement 
C4, 0.02-in-thick reinforcement 




Fig. 7 Effects of cutout-reinforcement thickness on the radial-displacement response of a geometrically perfect 
compression-loaded cylindrical shell with a cutout. The reinforcement is a 2.4-in. square-shaped reinforce- 
ment with 0°-ply orientation. The location of the axial and circumferential traces are indicated in Fig 6a 



Fig. 8 Effects of cutout-reinforcement thickness on the predicted initial-buckling deformations of a geometri- 
cally perfect compression-loaded cylindrical shell with a cutout. The reinforcement is a 2.4-in. square-shaped 
reinforcement with 90°-ply orientation (The outer boundary of the reinforcement region is marked in the 
figures). 
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Cl, no reinforcement 
C5, 0.005-in-thick reinforcement 
C6, 0.01 -in-thick reinforcement 
C7, 0.02-in-thick reinforcement 




Fig. 9 Effects of cutout-reinforcement thickness on the radial-displacement response of a geometrically perfect 
compression-loaded cylindrical shell with a cutout. The reinforcement is a 2.4-in. square-shaped reinforce- 
ment with 90°-ply orientation. The location of the axial and circumferential traces are indicated in Fig 6a. 



Fig. 10 Effects of cutout-reinforcement thickness and reinforcement fiber orientation on the predicted load- 
shortening response of a geometrically perfect compression-loaded cylindrical shell with a cutout. (4.4-in 
square-shaped reinforcement) 
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-0.47 


w/t 


0.35 


0.57 


a) Shell C8, 0.005 -in-thick 


reinforcement 


-0.25 w/t 0.46 

b) Shell C9, 0.01 -in-thick reinforcement 


-0.40 w/t 

c) Shell CIO, 0.02-in- thick reinforcement 


Fig. 11 Effects of cutout-reinforcement thickness on the predicted initial-buckling deformations of a geometri- 
cally perfect compression-loaded cylindrical shell with a cutout. The reinforcement is a 4.4-in. square-shaped 
reinforcement with 0°-ply orientation (the outer boundary of reinforcement region is marked in the figures). 


Cl, no reinforcement 
C8, 0.005-in-thick reinforcement 
C9, 0.01 -in-thick reinforcement 
CIO, 0.02-in-thick reinforcement 




Fig. 12 Effects of cutout-reinforcement thickness on the radial-displacement response of a geometrically per- 
fect compression-loaded cylindrical shell with a cutout. The reinforcement consists of a 4.4-in. square-shaped 
reinforcement with 0°-ply orientation. The location of the axial and circumferential traces are indicated in Fiu 
6a. s 
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-0.67 


w/t 


0.39 


-0.45 


w/t 


0.53 


-0.37 w/t 0.45 


a) Shell Cl 1 , 0.005-in-thick reinforcement b) Shell C12, 0.01-in-thick reinforcement c) Shell C13, 0.02-in-thick reinforcement 


Fig. 13 Effects of cutout-reinforcement thickness on the predicted initial-buckling deformations of a geometri- 
cally perfect compression-loaded cylindrical shell with a cutout. The reinforcement is a 4.4-in. square-shaped 
reinforcement with 90°-pIy orientation (the outer boundary of reinforcement region is marked in the figures). 

Cl, no reinforcement 
C11, 0.005-in-thick reinforcement 
C12, 0.01 -in-thick reinforcement 
C13, 0.02-in-thick reinforcement 




Fig. 14 Effects of cutout-reinforcement thickness on the characteristic radial-displacement response of a geo- 
metrically perfect compression-loaded cylindrical shell with a cutout. The reinforcement consists of a 4.4-in. 
square-shaped reinforcement with 90°-ply orientation. The location of the axial and circumferential traces are 
indicated in Fig 6a. 
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Fig. 15 Effects of cutout-reinforcement thickness and reinforcement fiber orientation on the predicted load- 
shortening response of a geometrically perfect compression-loaded cylindrical shell with a cutout. (8.0-in 
square-shaped reinforcement) 



a) Shell C14, 0.005-in-thick reinforcement b) Shell Cl 5, 0.01 -in-thick reinforcement c) Shell C16, 0.02-in-thick reinforcement 


Fig. 16 Effects of cutout-reinforcement thickness on the predicted initial-buckling deformations of a geometri- 
cally perfect compression-loaded cylindrical shell with a cutout. The reinforcement is a 8.0-in. square-shaped 
reinforcement with 0°-pIy orientation (the outer boundary of reinforcement region is marked in the figures). 
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Cl, no reinforcement 
Cl 4, 0.005-in-thick reinforcement 
C15, 0.01-in-thick reinforcement 
Cl 6, 0.02-in-thick reinforcement 




Fig. 17 Effects of cutout-reinforcement thickness on the radial-displacement response of a geometrically per- 
fect compression-loaded cylindrical shell with a cutout, The reinforcement consists of a 8.0-in. square-shaped 
reinforcement with 0°-ply orientation, The location of the axial and circumferential traces are indicated in Fig 



a) Shell C 1 7 , 0 .005-in-thick reinforcement b) Shell C 1 8 , 0 .0 1 -in-thick reinforcement c) Shell C 1 9 , 0 .02-in-thick reinforcement 


Fig. 18 Effects of cutout-reinforcement thickness on the predicted initial-buckling deformations of a geometri- 
cally perfect compression-loaded cylindrical shell with a cutout. The reinforcement is a 8.0-in. square-shaped 
reinforcement with 90°-ply orientation (the outer boundary of reinforcement region is marked in the figures). 
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Cl, no reinforcement 
C5, 0.005-in-thick reinforcement 
C6, 0.01 -in-thick reinforcement 
C7, 0.02-in-thick reinforcement 




Fig. 19 Effects of cutout-reinforcement thickness on the radial-displacement response of a geometrically per- 
fect compression-loaded cylindrical shell with a cutout. The reinforcement consists of a 8.0-in. square-shaped 
reinforcement with 90 -ply orientation. The location of the axial and circumferential traces are indicated in 
Fig 6a. 
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